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Existence of global solutions for a Keller-Segel-fluid equations 

with nonlinear diffusion 
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Abstract: We consider a coupled system consisting of the Navier-Stokes equations and a 
porous medium type of Keller-Segel system that model the motion of swimming bacteria liv¬ 
ing in fluid and consuming oxygen. We establish the global-in-time existence of weak solutions 
for the Cauchy problem of the system in dimension three. In addition, if the Stokes system, 
instead Navier-Stokes system, is considered for the fluid equation, we prove that bounded weak 
solutions exist globally in time. 
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1 Introduction 


We study a mathematical model describing the dynamics of oxygen, swimming bacteria, and 
viscous incompressible fluids in M 3 . More precisely, we consider the Cauchy problem for the 
coupled Keller-Segel-Navier-Stokes system in R 3 x [0, T) with 0 < T 


(KS-NS) 


' dtn + u ■ S7n = A n 1+a — V • (y(c)reVc) , 

< dtc + u ■ Vc = Ac — At(c)n, 

„ dtu + t{u ■ V)re + Vp = Au — nV(j), divre = 0, 


( 1 . 1 ) 


where re, c, u and p are the cell density, oxygen concentration, velocity field and pressure of the 
fluid. Here a > 0 is a positive constant such that porous medium type equation of re is under 
our consideration. 

It is known that the above system models the motion of swimming bacteria, so called 
Bacillus subtilis, which live in fluid and consume oxygen. The functions y : M —> M and 
/t : R —> R represent the chemotactic sensitivity and consumption rate of oxygen. The constant 
r in the third equation is 0 or 1. When r = 1, u becomes the velocity vector of fluid solving 
Navier-Stokes equation. If the fluid motion is so slow, one may assume r = 0 so that u is the 
velocity vector satisfying Stokes system. 

The above system (11.11) was proposed by Tuval et al. in [T3] for the case a = 0, which 
can be extended to the case a > 0 when the diffusion of bacteria is viewed like movement in 
a porus medium (see e.g. El, 0, EU, HU, EE|). Throughout this paper, we call the above 
system a Keller-Segel-Navier-Stokes equations (KSNS) if r = 1 and Keller-Segel-Stokes system 
(KSS) in case that r = 0. Compared to the Keller-Segel model (KS) of porus medium type: 


(KS) 


dtn = Are 1+ “ — V • (yreVc), 
rdtc = Ac — c + re, 
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where x is a positive constant and r = 0 or 1, we emphasize that the chemical substance 
(oxygen) in (11.111 is consumed, rather than produced by the bacteria. The existence of global- 
in-time bounded weak solution for Cauchy problem to (KS) with arbitrary large initial data is 
guaranteed when a > 1/3 for three dimensional case. Otherwise, the solution may blow up in 
finite time unless sufficiently small initial condition is assumed (see the results [10] and [?]). 

The main concern of this paper is to specify the values of a so that the global-in-time weak 
solutions for (KSNS) and bounded weak weak solutions for (KSS) are established, under the 
some conditions of x and The notions of weak and bounded weak solutions mentioned above 
are defined as follows: 

Definition 1 (Weak solutions) Let a > 0 and 0 < T < oo. A triple ( n,c,u ) is said to be a 
weak solution of the system m if the followings are satisfied: 

(i) n and c are non-negative functions and u is a vector function defined in R 3 x (0, T) such 
that 

n( 1 + \x\ + |logn|) £ L°°(0, T; L 1 (R 3 )), Vn^ € L 2 (0, T; L 2 (R 3 )), 
n £ L°°(0,T;L p (R 3 )), Vn 2 ? G L 2 (0,T;L 2 (1 3 )), 1 < p < 1 + a 


cGL oo (0,T;i7 1 (R 3 ))nL 2 (0,T;i? 2 (M 3 )), c £ L°°(R 3 x [0,T)), 
u £ L°°(0,T;L 2 (R 3 )), Vu £ L 2 (0, T; L 2 (R 3 )), 

(ii) (ro, c, u) satisfies the equation «np in the sense of distributions, namely, 

/ / i—mpt + Vn 1+a • V<p — nu ■ V<p — nx(c)Vc • V<p) dxdt = / no<p(-, 0 )dx, 

Jo Jr 3 Jr 3 

/ / (—apt + Vc • V<p — cu ■ V<p + nn(c)(p) dxdt = / cop>(-,0)dx, 

Jo Jr 3 Jr 3 

(—u ■ ift + Vrt • V-0 + (t(u ■ V)u) ■ if + riS7(f> ■ if) dxdt = uq ■ if(-, 0 )dx 

Jr 3 


n 

to Jr 3 


for all test functions ip 6 Cq° (R 3 x [0, T)) and if 6 Cq° (R 3 x [0,T),R 3 ) with V • if = 0. 

Definition 2 (Bounded weak solutions) Let a > 0 and 0 < T < oo. A triple (n,c,u) is said 
to be a bounded weak solution of the system hl.l\ ) if (n, c, u) is a weak solution in Definition 
[7] and furthermore satisfies the following: For any p £ [1, oo) and q £ [2, oo) 

(i) n £ L°°(( 0, T) x R 3 ), Vn 2 ? € L 2 (0,T;L 2 (R 3 )). 

(ii) c£ Li(0,T-W 2 ’ q (R 3 )), c* £ L«(0,T;L«(R 3 )). 

(hi) u £ L\ 0, T; fT 2 ^(R 3 )) ; u t £ T 9 (0, T; L"(R 3 )). 

Before we state our main results, we recall some known results in case that a > 0 (compare 
to e.g. m , ®, a, m and m for the case that a = 0 and references therein). It was proved 
in [5l that weak solutions of (KSS) for bounded domains exist in case that a £ (1/2,1] in two 
dimensions or in case that a £ ( ~ 5+ 1 '^ 23 7 ,1] in three dimensions. It was also shown that if 
domain is the whole space, i.e. R 2 or R 3 , then weak solutions exists globally when a = 1. 
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In [8], the exponent a is reduced up to 1/3 for the case that spatial domain is R' 3 , under the 
following assumptions on x and k: 


x(c),k(c),x'(c),k'(c) > 0, k(0) = 0, 


(. x(c)k{c))' 

X(c) - ’ 



( 1 . 2 ) 


For the case of bounded domains in dimension two, bounded weak solutions of (KSS) are 
constructed in m for any a > 0. In case that fluid equation is the Navier-Stokes equations, 
it was proved recently in [3j Theorem 1.8] that bounded weak solutions of (KS-NS) exist for 
any a > 0 in R 2 . 

On the other hand, in three dimensional case, p]2j considered a special case of y(c) = 1 
and k(c) = c and showed that global weak solutions of (KSS) exist whenever a > 1/7 for 
bounded domains. For the case of R 3 , [3] proved that global bounded weak solutions of (KSS) 
exist when a > 1/4 under the hypothesis either x'(‘) > Xo > 0 or k'(-) > ko > 0, where xo 
and K0 are positive constants. It was shown, very recently, in m that in case that a > 1/6, 
bounded weak solutions are constructed in a bounded convex domain Q C R 3 with smooth 
boundary with more relaxed conditions on x and k in general setting (see [16, Theorem 1.1] 
for the details). 

The main purpose of this paper is to establish existence of global weak and global bounded 
weak solutions for the Cauchy problem of the system (KS-NS) in R 3 under rather relaxed 
conditions of x and k compared to (jl.2jl . and smaller value of a ever known. 

Before stating our results, we address some conditions for k. To preserve the non-negativity 
of the density of bacteria n(x,t ) and the oxygen c(x,t ) for 0 < t < T, we need the condition 
k( 0) = 0 (see the proof of Lemma[I|. The condition k(-) > 0 is also essential since the bacteria 
is consuming the oxygen. Namely, we assume the following hypothesis: 


«;(•) > 0 and k(0) = 0. 


(1.3) 


We also present two different types of further assumptions on chemotactic sensitivity x an d 
consumption rate k together with the range of a. The first one is related to weak solutions. 

Assumption 1 We suppose that x, k satisfy x' G L™ c and k € L^ c with (11.31) . We assume 
further that one of the following holds: 

(i) a > 1/6. 

(ii) a > 0 and x'(') — Xo for some constant xo > 0. 

(in) a > 0 and «'(•) > Ko f or some constant no > 0. 

Another hypothesis is prepared for bounded weak solutions. 

Assumption 2 We suppose that x, k satisfy x' £ and k € L(f )c with (11.31) . We assume 
further that one of the following holds: 

(i) a > 1/6. 

(ii) a > 1/8 and x'(') > Xo f°i" some constant xo > 0. 

(in) a > 1/8 and k'( ) > kq for some constant kq > 0. 
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We are now in a position to state the main results. The first main result of this paper is 
the existence of global-in-time weak solution of the system (KS-NS) with r = 1, which means 
that the fluid equations are the Navier-Stokes equations. More precisely, the first result reads 
as follows: 


Theorem 1 Let r = 1 and the Assumption \7\ hold. Suppose that initial data (no,co,uo) 
satisfies 

n 0 (l + |x|-|-| log n 0 |) 6 L 1 (M 3 ), n 0 G L 1+ "(I 3 ), c 0 € L°°(M 3 ) n i? 1 (M 3 ), u 0 G L 2 (M 3 ). (1.4) 


Then, for each T > 0, there exists a weak solution (n, c , u) for the system ed and it satisfies 


sup 

o <t<T 


n{t) (|log n(t)\ + 2(x)) dx + ||n(t)||^t+ Q + ||Vc(t)||| 2 + \\u{t )\\ 2 L 2 


+ 


a 


l + Q 


Vn 2 (t) 


L 2 


+ 


1 _i_2 a 2 n r 

Vn— it) ^ + \\Ac(t)\\h + l|V«(t)||i 2 )dt<C, 


where (x) = (l+\x\ 2 ) 2 andC = C(T , ||n 0 || L i nL i+ a , ||n 0 logn 0 || L i , ||n 0 (a:)|| L i , ||Vc 0 || L2 


(1.5) 

IKM- 


If the fluid equations are restricted to be the Stokes equations, i.e. r = 0, and if the range 
of a in the hypothesis in Theorem [T] is a bit more restrictive, we can construct bounded weak 
solutions for the Cauchy problem. More precisely, our second result reads as follows: 


Theorem 2 Let r = 0 and the Assumption^ hold. Suppose that initial data (no,co,no) satisfy 
(11.41) and for any q G [2, oo) 

n 0 € L°°(M 3 ), c 0 € u 0 G W 1>9 (R 3 ). (1.6) 


Then, for each T > 0, there exists a bounded weak solution ( n,c,u) for the system (ED such 
that it satisfies 


\n\ 


L°°((0,T)xI 


+ 


P+Q= 

Vn 2 


L 2 ((0,T)xK 3 ) 


+ c 


Li (0,T;W 2 ’ q 


+ ll^t c llL'J(o,T;L'J(R 3 )) + ll U llL'J(0,T;iy 2 ^(R 3 )) + \\dt U W Li(0,T;L<i(R 3 )) < 

where C = C(T, \\n 0 \\ L oo , ||co|| w i, g , \\u 0 \\ w i, q ). 


(1.7) 


Remark 1 The results in Theorem [H can be rephrased as follows: In case that a > 1 /6 ; the 
bounded weak solutions can be constructed with the assumptions that y, k satisfy x' £ 
k G L^ c with «(■) > 0 and k(0) = 0 (compare to Uffi for bounded domain case). If we assume 
further that either x'(‘) ^ Xo > 0 or n'(-) > hq > 0, then bounded weak solutions exist for 
a > 1/8, weaker than a > 1/6. 


Our main concern is construction of weak and bounded weak solutions for the Cauchy 
problems of the system (EH) in K 3 but it can be extended without so much difficulty to bounded 
domains with Neumann boundary conditions for n and c and no-slip boundary conditions for u. 
To be more precise, let li C K 3 be a bounded domain with smooth boundary and we consider 
the system (11.11) in H x [0, T) with boundary conditions 

= 0, u = 0 on On. (1.8) 

ou on 

We then obtain following results for the case of bounded domains by following almost same 
arguments as in Theorem [T] and Theorem [2j 
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Theorem 3 Let r = 1 and the Assumption\7\hold with replacement of R 3 by a bounded domain 
ft with smooth boundary. Suppose that initial data (no,co,uo) satisfies 

n 0 £ T 1 (fi) n L 1+a (n ), Co € L°°(fi) n G L 2 (n). (1.9) 

Then, for each T > 0, there exists a weak solution (n, c, u ) /or f/ie system (O) and (ESD. and 
it satisfies 


sup ff n(t)\logn(t)\dx + \\n(t)\\]^+ a + \\Vc(t )\\ 2 L 2 + \\u(t )\\ 2 L 2 

<t<T \J R3 


0<t<T 
rT 


+ 


>0 


l + a 


Vn 2 (t) 


L 2 


+ 


Vn^(t)| L2 + ||Ac(t)||2 2 + ||V«(t)||| 2 )dt<C, 


where C = C(T, ||n 0 || L i nL i+a , ||n 0 logn 0 || L i , ||Vco || L 2 , ||w 0 || L 2 ). 


Theorem 4 Let r = 0 and the Assumption hold with replacement of R 3 by a bounded 
domain ft with smooth boundary. Suppose that initial data (no,co,uo) satisfy (11.41) and for 
any q £ [2, oo) 

n 0 G L°°(Q), c 0 £ W 1 ’ 9 ^), u 0 £ W 1 ’^ ft). 

Then, for each T > 0, there exists a bounded weak solution (n, c, u) for the system HMD and 
(USD such that it satisfies 


n 


L°°((0,T)xfi) 


+ 


P + OL 

Vn 2 


L 2 ((0,T)xn) 




+ II^ c IIl'?(0,T;I,9(Q)) + \\ u \\li(0,T-,W 2 :1(CI)) + ll^llL^TjL^fi)) < C, 
where C = C(T , \\n 0 \\ LOO , ||c 0 || w i, g , HwqII^i,,). 


Since verifications of Theorem [3] and Theorem [4] are similar to those of Theorem Q] and 
Theorem O only differences are indicated in the proofs in section 4. 


Remark 2 Theorem [7] and Theorem^ are improvements over the results in |3, Theorem 1.5 
and Theorem 1.7]. We also remark that Theorem [/] covers the result of ]12f . since the case 
X = 1 and k(c) = c is the special case of (in) in the Assumption d and besides, a > 1/8 is 
wider than a > 1/7 in fW/ . 

This paper is organized as follows. In section 2 and section 3 a priori estimates of weak and 
bounded weak solutions are established. In section 4 we present the proofs of main results. 


2 Weak solutions 

Throughout this section, we study the solutions of the approximate problem of (H3D given by 

d t n Q + u Q ■ Vn g = A (n e + g) 1+a - V • (x(c g )n g Vc g ), 

< d t c g + u g ■ Vc e = A c e - K(c e )n g , (2.1) 

. d t u e + t{u q ■ V)u g + Vp g = A u g - n g V(j), div u e = 0, 
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in R 3 x (0, T) with smooth initial data (no g , cq q , uo g ) given by 


n 0g = ^e* n o, °0e = fi e *c 0 and u 0g = fi g *u 0 


where <f) g denotes the usual modifier with g £ (0,1). 

It is known that, due to the standard theory of existence and regularity as done in [5] and 
[12] . there exists a classical solution of the equation (12.11) locally in time for each g £ (0,1). 
The main objective of this section is to derive appropriate uniform estimates, independent of 
g, of the solutions. The estimates are crucially used in Section 3 to extend the above local 
solution to any given time interval (0,T) and to construct the weak solutions and the bounded 
weak solutions of the equation (11.11) . 

We start with some notations. For 1 < q < oo, we denote by IT fc,9 (R 3 ) the usual Sobolev 
spaces, namely VF fc,9 (]R 3 ) = {/ £ L 9 (M 3 ) : D a f £ L 9 (R 3 ),0 < |a| < k}. The set of q— th power 
Lebesgue integrable functions on R 3 is denoted by L 9 (R 3 ). In what follows, for simplicity, 
|| • ||p denotes || • || iP ( R 3) for 1 < p < oo, unless there is any confusion to be expected. We 
also denote by W~ k,q (R 3 ) dual space of Wq ’ 9 (R 3 ), where q and q' are Holder conjugates. The 
letters C = C(*, ..., *) and C' = C'{*, ..., *) are used to represent generic constants, depending 
on which may change from line to line. 

We first recall that maximal regularity estimate of the inhomogeneous heat equation, which 
we use later. Let 0 < T < oo and 1 < p < oo, and we consider 

vt — Au = / in R 3 x (0, T) 

with initial datum v(x,0) = vo(x) with vq £ kF 1,p (R 3 ). Then following LP estimate is well- 
known: 


ll' u i|lz,P((0,T)xR 3 ) + IMIlp((0,T);W 2 ’P(R 3 )) — ^ \ II/Hlp(( 0 ,T)xR 3 ) + II u 01 |w 1 -p(R 3 ) ) ■ (2-2) 

In the following lemma, we give an estimate of solutions of (12.11) under the Assumption Q] 
The estimate is used in Section 3 to construct the weak solution of the equation (ED. For 
the sake of simplicity, throughout Section 2, we denote n g , c g and u e by n, c and u. Also, we 
define the functionals EM(t ) and D{t) as follows: 


E M (t) ■■= f n(t) (|logn(f)| + 2(x}) dx + ||n(t)||^“ + ||Vc(t)||| + 
J R 3 


M + 2 




and 


1 +oc 


D(t ) := Vn 2 (t) 


+ 


Vn^(t )\\ 2 + \\Ac(t)\\l + \\Vuml 


where (x) = (1 + |x | 2 )2 and M is a positive constant, which will be specified later. 


(2.3) 

(2.4) 


Lemma 1 Let T > 0. Suppose that (n, c, u) is a classical solution for the system \2.1\) c . 
g £ (0,1) with the smooth initial datum (no g , co Q , uq q ) satisfies the initial condition jl-4\) in¬ 
dependently of p. Assume further that x, k and a satisfy Assumption [7J Then, there exists 
C > 0 and M > 0, which are independent of p, such that for any 0 < t < T, EM{t) and D(t ) 
defined in (12.31) - (12.41) satisfy 


sup Em{t) + f D(r)dT < C, 
0<T<t J 0 


(2.5) 


where C = C(T , ||n 0 || L i , ||n 0 logre 0 || L i , ||n 0 (x)|| L i , ||n 0 || L i +a , ||Vc 0 || L 2 , IKIIl 2 )- 
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Proof. We observe, by integrating both sides of (12. lib . that the total mass of n is pre¬ 
served, i.e. ||ri(t)|| x = 11TT-o11! • It is also obvious by applying maximal principle to (12.lib that 
||c|| LO o (R 3) < ||coHoc j where M;/ := M 3 x [0,T). We also note here that n and c preserves non¬ 
negativity of the initial data by the condition n(0) = 0 and the parabolic comparison principle. 

• Case (i) in the Assumption [l] Here we consider the case that a > 1/6. We separate 
the range of a into three parts, that is, 1/6 < a < 1/3, 1/3 < a < 1 and a > 1, and treat each 
case individually. We start with the case that l/6<a<l/3. 


(Case 1/6 < a < 1/3): Multiplying (12.lib with (1 + logn) and integrating it parts, 

/- [ n log note + f V log n • V(n + g) 1+a dx = / Vn • (x(c)Vc) dx. 
dt J r3 Jr3 Jr3 

Since V log n ■ Vn = 4 | Vn 1//2 | 2 , we have 

/ Vlogn • V(n + g) 1+a dx = / Vlogn • (1 + a)(n + g) a X7n dx 


> / Vlogn • (1 + a)n“Vn dx = 


1 -|- (X 


1 +a 


Vn 2 


Taking into account 1/6 < a < 1/3, which admits (1 — a)/2 > 0, we obtain 


if 

dt J re 


< 


n log ndx + 

r 

1 + a Jm 


1 + a 

1 + Q- 

Vn 2 


1 + a 


Vn 2 


< / Vn • (x(c)Vc) dx 


1 — a . 

n^~ | Vc|) dx, 


where x := su Pk|, \x ( c ('))I • Applying Young’s inequality, we observe 


1+a 

Vn 2 

(n 2 Vcl ) dx < e / 

1+a 

Vn 2 


^ ’ Jr 3 



dx + C(e) / n 1 “ |Vc| dx. 


Combining above estimate and choosing sufficiently small e > 0, we have 


d_ 

dt 


n log ndx + C' 


1+a 


Vn 2 


< C I n 1_ “ |Vc| z dx. 


( 2 . 6 ) 


Reminding that 1 — 3a > 0 (since a < 1/3) and that c € L°°(Mj+ the term f R3 n 1 a | Vc| 2 dx 
is estimated as 


r Vc ■ Vc dx 


[ n 1 “ |Vc| 2 dx = [ n 1 
Jr 3 Jr 3 

C[f |Vn 1_ “| |Vc| dx-\- f n 1_ " |Ac| dx) 

\J r 3 Jr 3 J 


< 


< 


1+q 


Vn 2 


+ 


C(e)n 1 ~ 3a |Vc| 2 )dx + C [ n 1_a |Ac| dx. 

Jr 3 


Choosing small e > 0 in the above, we estimate (12.611 as follows: 


d_ 

dt 


[ nlogn dx + C\ Vn 1 ^ 2 < C x [ n 1 " 30 |Vc| 2 + In 1 '"! |Ac| dx (2.7) 
Jr s 2 Jr 3 
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for some C\ > 0 and C[ > 0. Next, Multiplying mh with n a and integrating it by part gives 


1 d 1+a 4a(l + a ) 


_||77<|| ' —I— _i_ 

1 +adt n 1+a (l + 2a) 


1+2 a: 


Vn 2 


= / Vn Q • nx(c)Vc dx 


< 


2«X 


„ l + 2a 

Vn 2 


n 2 IVcl | dx 


1 + a j_ 

Applying Young’s inequality to the integrand in the right-hand side of the above, we have 


d 


11+a 


J t \\ n \\i+a + C ' 


l+2c* 


Vn 2 


< C 


[ n |Vc| s 

J R 3 


dx 


( 2 . 8 ) 


for some C' > 0. Since the integral term of the right-hand side in (| 2 . 8 I) is estimated as 
/ n|Vc | 2 dx= nVc-\7cdx 

J r 3 Jr 3 

<C[ |Vn| |Vc| dx+ n|Ac| dx 


< C 


n 2 


l + 2a 


Vn 2 


|Vc| dx+ n|Ac| dx 


it follows from Young’s inequality that 


jL ||r?ll 1+a 4- r' 

dt ll n lll+« + C 2 


l+2a 


Vn 2 


< 


C 2 (J n 1-2a |Vc | 2 + n |Ac| dx\ (2.9) 


for some C 2 > 0 and C 2 > 0 . On the other hand, multiplying ( 12 . 11) 9 with —Ac and using 
integration by part, we get 


ld_ 
2 dt 


|| Vc ||2 + || Ac ||2 < / (n • Vc) Ac dx + k / n|Ac| dx 
Jr 3 Jr 3 

where k := max t 3 |n(c(-))| . Since the term f R3 (u ■ Vc) Ac dx in the above is estimated as 

(« • Vc) Ac dx = T / UiC Xi c x . x .dx = - V] / u ix .c Xi c x . dx 
l<t+3 y « 3 1£I+3' /k3 

^ ^ ^i,xjCC X j Xi dx C ||Vn || 2 ||v c11 2 


/ 


i<+i<3 ‘ 


where rq and c a , i denote the i-tli component of u and respectively, we have 


d 

dt 


IIVc||| + ||Ac||l < C 3 f ||Vu || 2 || Ac|| 2 + f n |Ac| dx 


( 2 . 10 ) 


for some C 3 > 0. For the fluid equation, energy estimate gives 

— ||n || 2 + 2 || Vn || 2 < C 4 / n|u| dx (2-11) 

dt J R 3 

for some C 4 > 0. Finally, to bound f n |logn| in d23D, we recall that (see e.g. [ 8 ] and El) 

/ n|logn|dx< / nlogndx + 2 / ( x)ndx + C , (2-12) 
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where (x) = (1 + |x | 2 ) 1//2 and (see e.g. (18) of [3]) 

[ (x)ndx = f nu-V(x)+ [ (n + g) 1+a A(x) + f V(x)-nx(c)Vc 

dt J R3 J R 3 J R3 J K3 

<(7(1 + ||«|| 2 + 11Vc|||) + ^C(e) + e||Vn 2 ||^ . 

Summing up (12.7IW(I2.13|) with sufficiently small e > 0, we have 

(/ R 3 n(l°gn + 2 (x)) dx+ j|n||{+“ + ||Vc||| + \\u\\l\ 


d_ 

dt 


(2.13) 


(2.14) 


+ CL 


/ 


2 

l+2a 


Vn 2 

+ 

2 

Vn 2 


+ || Ac||| + || V«||l 


<C 5 


l-3a 1 


n 


|Vc| 2 dx + f n 1 2Q |Vc| 2 dx+ f n 1 a \Ac\dx 

J R 3 il 3 

+ / n | Ac| dx + || Vn|| 2 || Ac|| 2 + / n\u\dx] . 

J R 3 J R 3 / 

:= C 5 (I + II + III + IV + V + VI). 

for some C' 5 > 0 and C 5 > 0. Due to 0 < 1 — 3a < 2/3 and 0 < 1 — 2a < 2/3, we note that 


I < 


and 


f R3 ^C(ei) + ein 3 ^ |Vc | 2 dx, if 1/6 < a < 1/3; 
||Vc|| 2 , if a = 1/3, 

/ 3 (C(e 2 ) + e 2 « 3 j |Vc | 2 dx, 


II < 


where the positive constants e\ and e 2 will be determined later. Hence, it follows from Holder 
inequality and Sobolev embedding that 


<) C^ei) ||Vc||| + ei Unollf ||Ac || 2 , if 1/6 < a < 1/3; 


l|Vc|| 


2 5 


if a = 1/3, 


(2.15) 


and 

II<C(e 2 )||Vc || 2 + e 2 ||n 0 |||||Ac|| 2 . (2.16) 

To estimate III, we apply Holder and Young inequality to have III < C(e) ||ZAc ||2 • 

Noting that | < < 2 (which is obtained from the condition 1/6 < a < 1/3), we obtain, 

from Young’s inequality and Gagliardo-Nierenberg inequality, the following : 


n 


i 2 —a 


=a<c>oiir 


1+4 a 
2+3a 


1 + Q! 

Vn 2 


6—6a 
2+3a 

2 


l+4a 

<C(e , )||n 0 || 1 2+3a + e 


1+a 

Vn 2 


2 

2 ’ 


Hence by choosing sufficient e' > 0, we have 


III < C(e 3 ) + e 3 


1 + q 

Vn 2 


2 

2 


+ e 3 ||Ac||2 


(2.17) 
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where £3 > 0 will be specified later. We estimate IV via a calculation similar to the above as 
follows. Holder and Young inequality give IV < C(e) ||n ||2 + e ||Z\c ||2 • Since | < 2 +ea < ^ v ^ a 
1/6 < a < 1/3, it follows from Young and Gagliardo-Nierenberg inequality that 


n 


l 2 <^ 


l+ 6 a; 

2+6o! 

1 


1+2 a 

Vn 2 


6 

2+6o: 

2 


l+6a; 

< C(e') ||n 0 |ir 3 “ +e' 


Vn 


l+2o! 

2 


2 

2 ’ 


Therefore, we have 


IV < C(e 4 ) + e 4 


Vn 


l+2o! 

2 


2 

2 


+ e 4 || A.c|||. 


(2.18) 


It is straightforward that V can be estimated as 


V<C(e 5 )||Vn||2 + e 5 ||Ac||*. 


(2.19) 


The positive constants e 4 in (12.18[i and £5 in ()2.19[) will also be specified later. Finally, we 
estimate VI as follows. Via Holder and Young inequality and Sobolev embedding, we obtain 

VI < C ||n|| 6 |M|g < C(e) ||n|| 6 + £ ||V«|||. 

Since 1 < | < 3+6a and 0 < 2 + 6a < 2, ||n|| 6 is estimated, via Gagliardo-Nierenberg inequality 
and Young’s inequality, as 

(5)(l+2, 

Ml < \\noh 2+6 ' 

5 

Thus, we conclude that 


l + 2a 

Vn 2 


2 + 6 a 
2 


< C(e') + £' 


l+2o; 

Vn 2 


VI < C(e 6 ) + £ 6 


1+2 o- 


Vn 2 


+ £6 ||ViaH 2 • 


( 2 . 20 ) 


Summing up (I2.15p ^ (|2.20l) with choosing sufficiently small £ 4 ,..., £6 > 0, we finally have 


d_ 

dt 


(^3 n( ' 1 ° gn + dX + + H VC II2 + IMI 2 ) 

+ Cq 


( 2 . 21 ) 


/ 

1 + Q 

2 

l+2o: 


Vn 2 

+ 

2 

Vn 2 


+ ll^clll + IIVtaIH 


<c 6 


2 + 11^112 + 1 


for some C' 6 > 0 and Cq > 0. To absorb the term C§ ||Vxt||§ in the right-hand side of the above, 
we test u to the ()2.1H ^ and then multiply sufficiently large constant M > Cq to have 


™j t \\u\\ 2 2 + M\\Vu\\ 2 2 <CM 


n \u\ 


dx. 


( 2 . 22 ) 


As ( 12 . 201 ) . the integral in the right-hand side of the above is estimated as 



n |u| dx < (7(e) + e 


Vn 


l+2o: 

2 


Ve||Vn|| 


2 

2 • 


10 































































Hence substituting it to (I2.22|) . choosing sufficiently small e > 0 and then adding it to (12.211) . 
we have 


d_ 

dt 


n(logn + 2(x)) dx + ||n||*+“ + ||Vc|| 2 + 

2 


+ C' 


2 + 2 




/ 

1+a; 

2 

1+2 a 


Vn 2 

+ 

2 

Vn 2 


+ l|Ac|| 2 + Il^ u ll2 


< c 


12 + l ) ' 

Integrating in time variables and combining estimate (12.121) . we deduce (12.51) . 

(Case 1/3 < a < 1): Performing similar calculations as in the above case, we obtain 

d 
dt 


d_ 

dt 


<C 2 , e 


n 


1—a 


| Vc| 2 + n |Ac| dx ] + e 


1 + 0 


Vn 2 


j t IIVc||a + ||Ac||a < C 3 (j|Vu|| 2 || Ac|| 2 + n |Ac| dx ^ , 
M| 2 + 2 ||Vu|| 2 < C 4 / n|n| dx 

JR 3 


d_ 

dt 


and 


f (x)ndx = [ nu-'V(x)+ f (n + g) 1+a A(x) + ( V(x)-nx(c)Vc 

dt J R3 JR3 Jr3 Jr3 

<(75(1 + ||u||i + ||Vc||i) + ^C(e) + e||Vn 2 || 2 ^ 

Summing up (I2.24p ^ (|2.28l) with sufficiently small e > 0, we have 

j t Qf g n(logn + 2(x)) dx + ||n||{+“ + ||Vc|| 2 + \\u\\fj 


/ 

1 + 0! 

2 

l+2a 


Vn 2 

+ 

2 

Vn 2 


+ ^6 y 

<Ce( f n 1_a |Vc| 2 dx + 
\J R3 

= : C 6 (I + II + III + IV). 


+ l|Ac|| 2 + ||Vu|| 


n 


|Ac| dx + ||Vn|| 2 ||Ac||, + / n\u\dx 

Jr 3 


(2.23) 


/ n log n 

dx + C( 

Vn 1 ? 2 < Ci [ n 1-a |Vc| 2 dx, 

(2.24) 

Jr 3 


2 JR3 

+ <4 

l+2o 

Vn 2 

2 <C [ n|Vc| 2 dx 

2 ./»3 

(2.25) 


(2.26) 

(2.27) 


(2.28) 


(2.29) 


for some C' 6 > 0 and C6 > 0. Since 1/3 < a < 1, it is direct that 0 < 1 — a < 2/3 and 
0 < 2 + 6 a < and tl lus ) we obtain the same estimations as those in the previous case for 
I~IV. To be more precise, we have 


I < 


C+i) IIVc||| + ei ||no||i ||Ac|||, if 1/3 < a < 1; 


l|Vc|| 2 , 


if a = 1, 


(2.30) 
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II < C( 62 ) + €2 


l + 2a 


Vn 2 


+ C 2 || Ac || 2 , 


m<C(e3)||Vu||2 + e 3 ||A C ||2 


and 


IV < C(e 4 ) + 64 


l+2g 


Vn 2 


+ €4 ||Vu ||2 • 


Following the same procedure as that in the previous case, we conclude that 


d_ 

dt 


( [ n (logn + 2 (x)) dx + ||n||}+“ + ||Vc||| + 
\J R3 

+ C' 


2 ]\d + 2 


\u\ 


(2.31) 

(2.32) 

(2.33) 

(2.34) 


f 

l+a 

2 

l+2a 


Vn 2 

+ 

2 

Vn 2 


+ l|Ac ||2 + || V?x||| 


< C 


+ 1 ), 


for some M, C, C' > 0, which gives the result after time integration together with estimate 

CM). 


(Case a > 1): As in the previous cases, we obtain 


-[ 

dt Jk 3 


n log ndx + C[ 


1 + q 


Vn 2 


< / Vn • (x(c)Vc) dx 


(2.35) 


d „ 

— \\ n 
dt 

d_ 

dt 


I l+a 1 /-v 
I l+a + °2 


l+2ct 


Vn 2 


<Ci / n|Vc| 2 + n|Ac| dx, 

7r3 

< C 2 [ n|Vc| 2 dx fsee (12.81) 1. (2.36) 

2 .7r3 

l|Vc|| 2 + || Ac||| < C 3 (j|Vu|| 2 || Ac|| 2 + J n |Ac| dx^j , (2.37) 

d p 

— \\uW 2 + 2 ||Vu ||2 < C 4 / n|u| dx (2.38) 

dt J R 3 


and 


d_ 

dt 


(x)ndx = / nn-V(x)+ / (n + £>) 1 + “A(x) + / V(x)-nx(c)Vc 
1 J R3 JR3 7r3 

< C 5 (l + ||n ||2 + ||Vc|| 2 ) + ^C(e) + e||Vr 


l+a 

'n 2 


(2.39) 


for some constants C+C^Ci ~ C+ > 0. Summing up (|2.35p ^ (l2.39l) with sufficiently small 
e > 0 , we have 


d_ 

dt 


<C 5 


/ n(logn + 2 (x)) dx + 
J R3 

+ C 5 


n 


I l+CK 
I l+a 


+ II Vc|| 2 + |M| 2 


(2.40) 


/ 

l+a 

2 

l+2a 


Vn 2 

+ 

2 

Vn 2 


+ II Ac||| + ||Vn || 2 


n 


|Vc| 2 dx+ [ n | Ac| dx + j| Vn || 2 || Ac|| ; 
• 7 r 3 


+ / n|u|dx 


:= C 5 (I + 11 + HI + IV) 
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for some C' 5 > 0 and C 5 > 0. We now estimate I as follows. Taking into account that ++ > 1 

1-1-a: 

via a > 1, we have n < C(e) + en 2 due to Young’s inequality and, therefore, we get 

I < C(e) ||Vc||| + e [ n^r- |Vc | 2 dx. 

Jr 3 


The second term in the right-hand side of the above is estimated by 

/ n^~Vc-\/c dx <C 
Jr 3 

<C 


1 + Q 1 + q 

Vn 2 -Vc + n 2 Ac dx 


1 + Q! 


Vn 2 


+ IV c ll2 + ll n lll+a + Il^ c ll2 


We also note, due to Gagliardo-Nierenberg inequality and Young’s inequality that 


2+2 a 

I II 1+Q! ^ s-i II 11 2-i-3o: 

Nll+a — ^ INIIl 


1 + Q 

Vn 2 

6q 

2+3 “ < c’ + 

1 + Q 

Vn 2 


2 



where we used that 2 +d, a < 2. Combining estimates, we conclude that 


I < C(e i) 


+ l +ei 


1 + Q! 


Vn 2 


+ l|Ac|| 


(2.41) 


We observe that the other terms II ~ IV are estimated exactly the same as those in previous 
case. More precisely, we have 


II + C (€ 2 ) + £2 


l + 2a 


Vn 2 


+ 62 II Ac|| 2 , 


m<C(e3)||Vu || 2 + e 3 ||Ac|| 2 , 


IV < C(e 4 ) + e 4 


1+2q 


Vn 2 


+ e 4 || Vn 


Following similar procedures as before, we conclude that 

( f nlogn dx + ||n||}t“ + ||Vc||| + 

\J R3 


d_ 

dt 


2 , M + 2 2 

-— uL 


(2.42) 

(2.43) 

(2.44) 

(2.45) 


+ c 


/ 

1 + Q 

2 

1+2q 


Vn 2 

+ 

2 

Vn 2 


+ II Ac||| + HVnll 2 


< C 


+ 1 


for some M, C, C' > 0, which implies (12.511 . Since its verifications are just repetition of pre¬ 
vious computations, the details are omitted. 


• Case (ii) and (iii) in the Assumption |T] Here, we only deal with only the case (iii), 
since the case (ii) can be proved in almost the same manner. Due to the previous result, since 
the case a > 1/6 is direct, it suffices to consider the case 0 < a < 1/6. Multiplying (12. IIP with 
—Ac and integrating it by part gives 

||Vc || 2 + ||Ac || 2 < [ (n • Vc) Ac dx + [ n(c)nAc dx. 

2 dt JR3 J R3 
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Since the last integral of the above is estimated as 


[ n(c)nAc dx = - / (V(c)nVc + /e(c)Vn) Vc dx 
J R 3 Jr 3 

< — kq / n|Vc| 2 +77 / |Vn| |Vc| dx, 


we have 


d_ 

dt 


||Vc || 2 + ||A.c||| + 2n 0 [ n |Vc | 2 dx 

JR 3 


< C I || Vu || 2 || Ac || 2 + 


1 + Q 

Vn 2 


n 


^ IVc|) 


dx 


Multiplying (I2.1l) i with (1 +logn) and integrating it by parts, we get 


d f , 4 

— In log ndx H- 

dt Jm 3 1 + a 


1+a 


Vn 2 


< 


Vn • (x(c)Vc) dx 


< 


2 x 

1 + a 


1 + a 

Vn 2 


1 —a 

n 2 Vc dx. 


(2.46) 


(2.47) 


Since the integral in the last inequality of (12.461) and (12.471) is estimated, due to Young’s 
inequlity, as 


/ 

Jr 3 


1 + a 

/ 1 —a A 

1+a 

2 f 

1 —a 

Vn 2 

( n 2 Vc J < e 

Vn 2 

+ C( e ) / 

2 JR3 

n 2 Vc 


we obtain by choosing a sufficiently small e 
n log ndx + ||Vc|| 


d_ 

dt 



l+a 

+ 

Vn 2 


< 


+ ||Ac || 2 + 2ko / n|Vc| 2 dx 


C ( || Vn || 2 || Ac || 2 + n 1_ " | Vc | 2 dx) . (2.48) 


Multiplying (12.11) i with n“ and integrating it by parts gives 


I i+a , 4a(l + «) 
1 + a " 1+a (1 + 2 a ) 2 


n 


l+2ct 


Vn 2 


= / Vn Q ■ nx(c)Vc dx 


< 


f 

1 + a Jr e 


_ l + 2a 

Vn 2 


n 2 IVc| | dx. 


Young’s inequality implies 


^ II ||i+® 

1 +adt 


Vn 


1+2 a: 
2 


2 

2 


<c 0 


n |Vc | 2 dx 


for some Co > 0. Now we choose sufficiently small eo > 0 satisfying eoCo < 2 ko- Multiplying it 
by the both sides of the above gives 


ep d 
1 + adt 


n 


1+a 

1+a 


+ ep 


Vn 


l+2a 

2 


2 

2 


< epCp 


n |Vc | 2 dx. 


(2.49) 


14 























































Summing up (12.481) and (12.491) . we have 


d_ 

dt 


n log n + || Vc ||2 + 


£q 

1 OL 


\n\ 


1 + 0 ! 

l-\-a 



l + o: 

2 

1+2 a 

+ 

Vn 2 

+ eo 

2 

Vn 2 


+ l|Ac|| 2 + *0 


<C ( ||Vu || 2 ||Ac || 2 + / n 1 ~ a |V<+ dx 


n |Vc | 2 dx 


where <5o := 2no — cqCq > 0. Therefore, applying Young’s inequality n 1 a < C(e) + en with 
sufficiently small e > 0 , we obtain 


d 

dt 


nlogndx + ||Vc||p H--— ||n 

1 + a 



1 +a 

2 

l+2a 

+ 

Vn 2 

+ e 0 

2 

Vn 2 


+ IIA4 + I 


<0 (||Vc||l+ IIV^Ill) . 


n |Vc | 2 


Finally, absorbing the term ||Vzt||a as we have done previously (see (12.221) 1 and combining 
estimate (12.131) . we deduce (12.51) . This completes the proof. □ 


3 Bounded weak solutions 

In this section, we provide some uniform estimates of solutions to the system ([2.11) with r = 0, 
which is used in Section 4 to construct the bounded weak solution of the equation (11.11) and 
prove Theorem [2j To be more precise, our result is read as follows: 

Lemma 2 Let r = 0 and T > 0 be given. Suppose that (n, c, u ) is a classical solution for the 
system ( E 3 P „ q £ ( 0 , 1 ) with the smooth initial datum (no e , cq s , uq q ) which satisfies the initial 
conditions (HI and m- Assume further that and a satisfies Assumption [H Then, for 
any 0 < t < T 

n G L°°(0, T; L P (M 3 )), 1 < p < oo, (3.1) 

Vn 2 ? G L 2 (0,T;L 2 (M 3 )), 1 < p < oo, (3.2) 

c,u G ^(O.TjW 1 ’ 9 ^ 3 )) nL q (0,T-,W 2 ’ q {R 3 )), 2 < q < oo, (3.3) 

Ct, Ut G L q (0, T; L 9 (R 3 )), 2 < q < oo. (3-4) 

Proof. Since Assumption [2] is stronger than Assumption [U it is automatic, due to Lemma HI 
that 

n € L°°(0,T; L P (R 3 )) and Vn 2 ? G L 2 (0,T;L 2 (1 3 )), 1 < p < 1 + a. (3.5) 

In what follows, we complete the proof by showing that for each a > 1/8, if n satisfies (13.51) 
then m (13.4p hold. We remark that the first case (i) in the Assumption [ 2 l i.e. a > 1/6, 
cannot be relaxed as a > 1/8, since (13.51) is not available in case only a > 1/8 is assumed. 
First, we consider L p estimate of (I2.1D i 

--T-NIS+ [ Vn p_1 • V(n + g) p dx = — [ n p_1 V ■ (x(c)nVc) dx. 

Pdt P ./ R 3 J R3 
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Since Vra p 1 • Vn = (4 (p — 1 )/p 2 ) |Vn p / 2 | 2 > 0, we have 
1 d 


P+OC 

2 

Vn 2 

< 

2 


_.. up | 4(p-l)(l + a) 

pdt p (p + a ) 2 


Since the righthand side of the above is estimated as 

I n p ~ l V • (x(c)nVc)dx < C f 
J R 3 Jr 


/ n p 1 X7 ■ (x(c)nVc) dx 

J K3 


<c[ 

p+a 

Vn 2 

Jr s 



[J — CJt 

n 2 \\/c\dx, 


it follows from Young’s inequality that 

1 d 




V+ot 


Vn 2 


<Cp [ n p ~ a |Vc | 2 dx. 
J R3 


(3.6) 


Using Holder inequality, Sobolev embedding and Young’s inequality, we estimate the integral 
in the righthand side of the above as 


[ n p ~ a \Vc\ 2 dx < ||n p -“||_ E _ |Vc | 2 
J K 3 P~ a a 


< Cp\\n\\ p a l|Ac || 2 6 p <Cp(— + 


a p — a 


2p+3a 


\P 


P 


"IIS HAc|| 2 _^. 


2p+3o: 


Combining estimates, we get 

d 


dt 


< <Cp 3 (z + L ^ INI?) IIMLse. 


P P 
2 


2p+3a 
2 II A _11 2 


< Cp A ||Ac|| 6 P ||n||p + Cp 2 IIAc 


2p+3o: 


6 p 

2p+3a 


Thus it follows from Gronwall inequality that 


|n|| p < exp ( cp 1 


/W)||^ f \\AcW\p 

J 0 2p+3a J Jo 2p+3a 


ds + ||n 0 || p 


2p+3o: J Jq 2p+3o: 

and hence (13.11) holds except for p = oo whenever the following holds: 

f T 

l|Ac(,)|| 2 _ 6 p ds < oo, 1 d - a < p < oo. 

Jo 2p+3a 

We derive (13.71) by treating two cases, i.e. a > 1/3 and 1/8 < a <1/3 separately. 

(Case a > 1/3): Applying maximal regularity estimate ( 12 . 21 ) of heat equation to ( 12 . 11) 9 

[ ||Ac(s)|| 2 6 P < C (1 1Vcp1 1 2 6 P + f ||n(s )|| 2 6 P ds + [ \\u • Vc || 2 e P ds 

JO 2p+3a \ 2p+3a J Q 2p+3a J Q 2p+3a 


■=C[ 11Vc 0 11 6 P +I + II). 

2p+3a 


(3.7) 
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Now, for p > 1 + a, we obtain via interpolation 


I < C 

< c 


l 

i: 


T 0 (l + 2g)(4p-3a) (l + 2a) (4p-3a) 

l"Wlir IlnWIlilS™ 1- * 


_ l+2a 

Vn 2 (sj 


4 _ 3a 

l + 3a p(l + 3a) 


ds = C 


f 


_ l+2a 

Vn 2 (s) 


2-S p 

ds 

2 


for some d p > 0. We note that the last equality of the above holds becausea > 1/3. Before 
estimating II, we first observe that 


u £ L°°( 0 , T ; L 6 (M 3 )). 


(3.8) 


Indeed, since r = 0, the vorticity, ui = V x u, satisfies ut — Aw = —V x (nV<^>), and its L 2 
estimation becomes 


— ||w ||2 + ||Voj|| 2 < C f Vw • ndx < e || Vow||^ + C(e) ||n || 2 
dt J R3 


Hence we have 


d_ 

dt 


k jj 


\1<C\ 


n 


o_ 3+3o: 

l2<GINIi^ +3a 


1 +a 

Vn 2 


2+3a 

2 


(3.9) 


Since 0 < 2+ 6 3a < 2 due to a > 1/3, we obtain f)3.8|> by Sobolev embedding. Next we 
estimate II. Applying maximal regularity estimate (12.21) for heat equation and (13.81) . we have 
for p > 1 + a, 


II < 


< 


r T r T 

/ IKs)ll <5 ll Vc ( s )ll!_ 6 ^ ds < C / ||Ac(s)|| 2 2 P ds 

JO P+3a Jq p+a 

c f||Vc 0 || 2 _ 2 p_ + / T (||n(s)|| 2 _^ + ||n • Vc(s)||\jd; 

\ P+a= Jq P+a p+a 


<C + C 


<c + c 


= :C + C 


o 3(l + 2a>)(p —a) 3(l+2a)(p —a) \ 

n(s )|| 1 2p(1+3a) ||ra( 8 )|| 3+ a ^ 1+3Ql) + ||Vc(s )|| 2 6p ) ds 

2p+3 a J 


3(l + 2a)(p — a) 


_ l + 2a 

Vn 2 (s) 


3 _ 3a 

l+3a p(l + 3a) 


+ ||Vc(s)|| 6 p ds 

2p+3a 


f 


l + 2a , s 2 — 5^ 


Vn—(s) + || Vc(s)||_ 6 p_ ds, 

2 2p+3a 


where 0 < d p : = + p ( 1 ^ 3 Q ,) < 2. We note here that for each q with 2 < q < 6 , 

f 0 T || Vc || 2 ds < oo due to the results in Lemma[TJ Therefore we conclude that n € L°°(0, T; L P (M 3 )) 

and VnV € L 2 (0, T; L 2 (R 3 )) for max{l + a, 3a} < p < oo, which can be easily extended 
tol + a<p<oo from (13.51) . It remains to show that n is bounded. Indeed, since 
n € L°°(0, V; L P (M 3 )) for all 1 < p < oo, we can see that c t , V 2 c, ut and V 2 n belong to 
L 9 ((0, T) x M 3 )) for all q < oo and therefore, we also note that Vc £ L°°((0, T) x M 3 )). Using 
the estimate (|3.6D and Vc £ L°°((0, T) x M 3 )), we obtain 


where (3 


d „ 
— N 
dt 


< Up 2 ||n||£_ a < Cp 2 ||n 


p-i 

l 


n 


p(p-a-l) 

I P-1 

IP 


< Cp 2 


InUf-W 


a/(p — 1). Via Gronwall’s inequality, we observe that 


n(t)|| p < (Cp 2 t)p + ||n 0 || p , t<T. 
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Passing p to the limit, we obtain for all t < T 

ll n ( i )llz,°°(]R3) <1 + ll^oll L°°(R 3 ) • (3.10) 

This completes the case that a > 1/3. 

(Case 1/8 < a < 1/3): We first show that 

n € L°°(0, T;L P (M 3 )) and Vn 2 ? G L 2 (0, T; L 2 (M 3 )), 1 < p < 1 + 4a (3.11) 

and we then derive (13.71) via (13.111) . To show (13.111) . we compute L p estimate of n as follows: 




'0 


Vn *2 (s) ds < f f n p 3 V ■ (x(c)nVc) dxds + 

2 Vo Vrs 


From Young’s inequality, the first term in the righthand side of the above is estimated as 

n n p_1 V • (x(c)nVc) dxds = — f f Vn p_1 • (x(c)nVc) dxds 

: 3 Vo Vrs 


Vn 2 


ds + C(ei) [ [ n p ~ a |Vc| 2 dxds. 
Vo Vrs 

We note that the last integral of the above is estimated as follows: 

[ [ n p ~ a |Vc| 2 dxds = f [ n p ~ a Vc • Vc dxds 
Vo Vrs Vo Vr3 

ft 


<C 


<C 2 


([ [ |Vn p "“||Vc| dxds+ [ [ |n p " Q | |Ac| dxds) 
Wo Vrs Vo Vrs J 


f 


P+CX 


Vn 2 


ds + C(e 2 ) 
2 



n p_3 “ |Vc| 2 dxds + C / [ \n p ~ a \ |Ac| dxds. 

Vo Vr3 


Hence choosing sufficiently small ei and 62 > 0, we have 


<C 


<C 


m p P + r 

Jo 

f f |tt p_3a | |Vc| 2 dxds + [ [ |n p_ “| |Ac| dxds ) + ||no|| p 
Vo Vrs Vo Vrs / 

ft 


3 “l \Vc\ 2 dxds + (||n|| p _“+} + ||Ac|| p _“+})ds) + 


:=C(I + II) + |K|| p . 

For each p with 1 + a < p < 1 + 4a, we estimate I via maximal regularity estimate (12.21) 


I < 


C f ||n p - 3 “(s)||^. |Vc(s)f 

JO p — 3a 


1 + q 


ds 


l+4a :—p 


<C j ||Ac(s)|| 2 i ds<C(||Vc 0 || 2 i + / (||n(s)|| 2 + ||u ■ Vc(s)|| 2 )ds ) , 


iri 


r 1 


irw 


(3.12) 
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where r± := • We note that Holder inequality is applicable to the first inequality of 

the above due to the conditions of a and p, that is, l/8<a<l/3 and 1 < p < 1 + 4 a. Now 
let us estimate the term fg ||u • Vc || 2 ds in (13.121) . Considering the following L 2 - estimation of 
the vorticity equation: 


4- 11 Cij 11 2 + HVwlla <c [ Vw • ndx < e \\Vuj \\ 2 2 + C(e) ||n || 2 , 

dt J m>3 


we have 


J „ 3(l+2cQ(l-cQ 

|lHl2<C'Nl2<C'IHI 1 ia 2+5Q 


1 + 2 Oi 

Vn 2 


3—6a 


2+5a 

2 


(3.13) 


Since 0 < < 2 due to 1/8 < a < 1/3, we have 


u e L°°(0,r ; L 6 (M 3 )). 


Taking into account that 1 < rq < 3 and the above, we note, due to Holder inequality, that 

/Vvcn>, riMisiivdi 6 r^ ds < C [ ||Ac|| 2 6 n ds 
Jo Jo e-n Jo e+n" 


+ ||u • Vc|| 6 ri )ds 

6+ri 6+r^ 


||Vc 0 || + [ (||n||V 

\ 6+r 1 J 0 6+r- 

<c( HVcoll v + [ (||n||+n_ + llVcII^Jds 

\ e+r 1 J 0 6+rj 


Substituting it to (|3.12l) and applying interpolation inequality, we obtain 


I < 


cfllVcollV +l|Vco|| r 2 i+ A||n|| 
V s+n Jo 


2(1~+) , I. I^C 1 —£>2) | U 1,202 

1+a Il' t ll3p+3a ~ r ll'dll 


"ll£U + ll v<)ds , 


where 9i = ^ P+ 4 ip+ 2 a- 1) +1 ^ an( I = 2 (i+a)( 3 fH- 3 cA) • We remark that it is not difficult to verify 
1 + cc < ri < 3p + 3a and 1 < < 3p + 3a whenever l/8<a<l/3 and 1 + a < p < l + 4a, 

which enables us to apply the interpolation inequality to the above. Therefore, we have 


(i 


1<C\ HVcollV +||Vc 0 || 2 + / ( 




6 + 7*1 


A 

p+a 

Vn 2 

49j 

p+a 

p+a 

Vn 2 

Jo 


2 



4^2 ^ 

2 P+ “ + II Veiled, 


=:C( HVcoll!^ + HVcoll^ + / ( 
V e + r i Jo 


p+a 

Vn 2 


5i 


+ 


P + Q 


Vn 2 


&2 


+ 11 Vc || 2 )ds 


Similarly, for each p with 1 + a < p < 1 + 4a, we estimate 
II < C [\\\n\Q + ||u||? HVcHl^M + ||Vcoll 

JO 6 -r 2 

<c( A||n ||? 2 + ||n||M + IIVeiled* + IIVcoIIM + IIVco|| 

\J o 6 +r 2 6 +r 2 


17*2 

Ir2 


r 2 

r 2 I ’ 


where r 2 := p — a + 1. We note that in the first inequality of the above, Holder inequality is 
applied due to the conditions of a and p, that is, 1/8 < a < 1/3 and 1 + a < p < l + 4a. Since 
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1 + a < r 2 < 3p + 3a, 1 < < 3p + 3a and 2 < r 2 < 6 , we obtain the following estimate 

via applying interpolation inequality: 

II < c(f\ ||n||g5 3 „ + ||n||? (1 - |,4) IMIg^J * 

+ /‘(IIVc||? ( 1 - 85) IIVeil? 9 ')* + HVcollV + IIVcolG ), 

J 0 6+^2 ' 

where 

_ 3(p + a)(r 2 - 1 - a) _ (p + a)(5r 2 - 6 ) _ 3(p - a - 1) 

3 " r 2 (3p + 2a —1) ’ 4 " r 2 ( 6 p + 6 a —2) &n 5 ' 2 r 2 

Therefore, we have 


II < C | 
= : C 


p+a 

Vn 2 


2i~2 fl 3 

p+a. 


+ 


p+q 

Vn 2 


2r2^4 

p+a 


+ IIAcie 8 -)* + lIVcolll^ + lIVcoll 


6+ r 2 


\r 2 

\V 2 


f 


P+a 


p+a 

Vn 2 

+ 

2 

Vn 2 


s 4 


||Ac||*)* + lIVcolll-, + lIVcoll 


6+1-2 


Since one can see that 0 < d* < 2 , i = 1 , 2 ,..., 5 whenever 1/8 < a < 1/3 and 1+a < p < l+4a, 
it follows from Young’s inequality that 


n(t)\\ P P + Ci 



p+q 

Vn 2 


2 

2 


ds < C 2 



ds + 1 


0 < t < T 


for p with 1 + a < p < 1 + 4a. We know that, for each q with 2 < q < 6 , Jo W^cfqds < oo. 

Since 2 < n < 6 whenever 2+110 < p < 1 + 4a, we have n e L°°(0, T ; L p ( M 3 )) and VnM € 
L 2 (0, T ; L 2 (M 3 )), for 2+310 < p < 1 + 4a, which can be easily extended tol + a<p<l + 4a 
from (|3.5D . Thus we conclude that (13.111) holds. 

We are now ready to derive (13.71) for the case 1/8 < a < 1/3. Let po = f Since 

1 < Po < 1 + 4a whenever 1/8 < a < 1/3, it is evident from the above conclusion that 


n e L°°(0,T;L Po (R 3 )). 


(3.14) 


From the maximal regularity estimate (12.21) for heat equation, we have 


/ T ||Ac|| 

Jo 


< c 


2p+3a 


|n|| _ 6 p ds + / ||n • Vc|| _ 6 p_ ds 


2p+3a 


2p+3a 


=: <7(1 + 11). 

6 p 

2p+3a 


For p > 1 + a, we have po < 2 J/+ < 3po + 3a and hence it follows from (13.141) that 


f'T o (Pn + a )(6P"2pp n -3qp n ) (PO+q)(6p-2pPQ-3qpo) 

I <c I l|n(s)|| Po 2(2po+3 “ ) ||n(s)||ZZI3^° T5 ^ ds 


' f T Ms 

Jo 

f 

f 


I 3po +3o: 


pp+q 

Vn 2 (s) 


12—4po 6app 

2p 0 +3a p(2pQ + 3a) 


ds 


PQ + <* . . 

Vn 2 (s) 


2 —<5r) 


ds 


20 




















































for some 5 P > 0. On the other hands, by the maximal regularity estimate (12.211 for the heat 
equation, we have for p > (in fact, p > 1 + a since 0 < < 1 + a) 


r T r T 

< / ||u(s)||g ||Vc(s) || 2 6 P ds<C / ||Ac(g)l| 2 2p ds 

Jo p+3q Jo p+“ 

< C f||Vc 0 ||^ 2 p_ + [ (||n(g)|| 2 2 P + \\u • Vc(s)|| 2 2 P )di 

\ P+Oi Jq p+ot P+ol 

( r T 2 —— 

<C ||Vc 0 || 2 ^+ / (||n(a)|| 1+a 

\ P+a Jo 

c( ||Vco|| 2 _^+ [ T ( 

V P+ a Jo 


2 _ 3(l+2a)(p— cap—ot — or) 


p(2+5a) 


K«)ll 


3(1+2a) (p —ap—a —a:) 


3+6a 


p(2+5a) 


< 


l + 2a 

Vn 2 (s) 


6—6 a: 6 a+6 a: > 

2 +5a p(2 +5a) + || Vc ( s )||2 jE _) ds 

2 2p+3a 


2 _ 

=: C ( ||Vc 0 ||j 2 p_ + / ( Vn^(s) P + l|Vc(s )|| 2 6p )ds 


P+ot 


2p+3 a 


+ II Vc(s)||^_6p_ 

2p+3a 



where 0 < 5' p := < 2. We remark that ||Vc|| 2 <is < oo for each q with 

2 < q < 6 , due to the results in Lemma [H 

So far, we have shown cud holds for each a > 1/8. Therefore we conclude that n € 
L°°(0, T; L P (R 3 )) for 1 + a < p < oo. Following similar procedures as in (13. 1011 , we can show 
that L°°-norm of n is bounded. Since arguments are on the same track, we omit the details, q 


4 Proofs of Theorems 

To prove Theorem [T] and Theorem [2j using the uniform estimates established previously, we 
construct weak and bounded weak solutions. Using the uniform estimates established in the 
previous section. Since the argument is rather standard (compare to PI 03 HU 121 HU), we 
omit the details and give the sketch of how our constructions are made instead. 


Proofs of Theorem [T] and Theorem [21 We consider only the case of Theorem [21 since 
the proof of Theorem Q] is essentially same. First, we recall the regularized system (11.111 with 
the initial data (no,, co,„ uo g ) which are chosen as smooth approximations of (no,co,uo): 


no e = * n 0 , c 0 , = ip e * c 0 and u 0g = i> g * u o, 


where 4> g denotes the usual molliher. The convergence of (no,, Co,, uo g ) entails that the es¬ 
timates obtained in Lemma Q] and Lemma [2] are uniform, independent of g, precisely, the 
constant C and M in (12.511 can be chosen independent of g. Likewise, there exists a constant 
C such that for q < oo 


Ml 


L°°((0,T)xI 


+ 


q+ot 

Vn , 2 


<C, 

(4.1) 

L 2 ((0,T)xR 3 ) 


i) + \\®t c Q Il9(0,T;L?(R 3 )) < 

(4.2) 

)) + \\dt u e L9(o,t ; L9(r 3 )) < C- 

(4.3) 


According to the estimates we have derived, a bootstrap argument can extend the local solution 
to any given time interval ( 0 , T) (compare to [5j, [E] and m for more detail). Let k be 
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any number with k > 2 + a. We then show that dtn Q and dtn k are, independent of g, in 
L 1 (0, T; W^ 2,2 (M 3 )), where IT“ 2,2 (M 3 ) is the dual space of IT 2,2 (R 3 ) (compare to [12]). Then 
via Aubin-Lions Lemma, by passing to the limit, we have some weak limit (n, c, n), which turns 
out to be a weak solution. Its verification is rather straightforward, and thus the details are 
skipped. It is also direct that (n, c, n) is a bounded weak solution and satisfies the estimates 
(14.1I) - (I4.3D . This completes the proof. □ 

Next we present the proofs of Theorem[3]and Theorem|4j As mentioned in the Introduction, 
we indicate only difference compared to the case of M 3 . 

Proofs of Theorem [3] and Theorem [4] We note that unlike R 3 , L 1 estimate of n(x) is 
not necessary, because negative part of J^nlogn is controlled by |MI_Li(m- We also observe 
that Gagliardo-Nierenberg inequality should be slightly modified, for example, we used in the 
case M 3 (see the inequality right above (12.181) 1 


n 


L 2 (R 3 ) > 


<C I 


n 


l+6a 
I 2 + 6 o: 


l+2a 

Vn 2 


2+6 a 

L 2 ( R 3 ) 


In the case of bounded domains, it is replaced by 


ra 


L 2 ( 0 ) ^ 


<C I 


n 


l+6a 
I 2+6a 

1x4(0) 


l + 2a 

Vn 2 


6 

2+6a 


L 2 (n) 


+ C\\n 


2 

L !( 0 ) • 


Major modifications lie in the estimate of vorticity, ca, because the boundary condition of uj 
is not prescribed. More precisely, in order to obtain (13.81) . we used the equation of vorticity, 
which is not useful to the case of bounded domains. Here we show (13.81) differently not by 
using vorticity equations. Let Q = x (0,T). Using L p —type estimate of the Stokes system 
(see e.g. El), we note that 

IKIIl 2 (Q) + \\ u \\l 2 ((0,T)-,H 2 {VL)) + \\^P\\l 2 (Q) ^ C (ll n ll L 2 (Q) + ll u oII m{o.)) ■ ( 4 - 4 ) 

Testing — An to the fluid equations, 

[ |Vn| 2 dx+ / |An| 2 dx<- f |A u\ 2 dx + C [ n 2 dx+ I \/S.u\\\7p\dx 
2 dt Jti Jq 4 Jq Jq 


< l i | Au| 2 dx + C [ n 2 dx + C [ \Vp\ 2 dx. 

Therefore, after integrating the above in time over (0,f) for any t < T and combining the 
estimate (|4.4I) . we obtain 


ll V “WllL2(n) + ll A-u llL 2 (Q t ) ^ II Vw o|Il2 ( Q ) + C (\\n\\ 2 L 2(Q t ) + ll V PllL 2 (Qi)) 

< C (jl'UollwA+n) + ll n lli 2 (Q t )) ) 

where Qt = (0,f) x ki and we used (14.41) . We note that L 2 norm of n can be estimated in the 
same ways as in (13.9 j) and (13.9(1 . which implies that Vn £ L°°((0, T); L 2 (fl)) and therefore, it 
is automatic that n £ L°°((0, T); L 6 (9)) via Sobolev embedding. The rest parts of proofs are 
essentially the same as the cases of whole space, and thus we omit the details. r-i 
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